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We incorporate single crystal FesCU thin films into a gated device structure and demonstrate the 
ability to control the Verwey transition with static electric fields. The Verwey transition temperature 
(TV) increases for both polarities of the electric field, indicating the effect is not driven by changes 
in carrier concentration. Energetics of induced electric polarization and/or strain within the FesC^ 
film provide a possible explanation for this behavior. Electric field control of the Verwey transition 
leads directly to a large magnetoelectric effect with coefficient of 585 pT m/V. 

PACS numbers: 75.47.Lx, 71.27.+a, 75.85.+t, 71.30.+h 



Electric field control of magnetic and metal-to- 
insulator transitions in highly correlated materials has 
generated great interest both scientifically and techno- 
logically [15 . Magnetite (Fe304) is a highly correlated 
material that undergoes the well-known Verwey transi- 
tion with sharp changes in the electric, magnetic, and 
structural properties [6 -8 . Following over seven decades 
of intense research since its discovery in 1939, the origin 
of the Verwey transition as a charge order-disorder tran- 
sition has become clear experimentally just in the last few 
years [9HT5]. In this Letter, we report the surprising ob- 
servation of electric field control of the Verwey transition 
in Fe304 thin films. An electric field applied by electro- 
static gates is found to stabilize the Verwey structure and 
increase the transition temperature (Ty). Furthermore, 
electrical control of Ty leads to a new mechanism for 
generating a magnetoelectric effect. We obtain a magne- 
toelectric coefficient of 585 pT m/V, which is one of the 
largest for a single-phase material. These results provide 
an alternative approach to advanced electronics, infor- 
mation processing, and spintronics [T6HT9] . 

Magnetite is the oldest known magnetic material (lode- 
stone) and the Verwey transition at Ty ~ 120 K is one 
of the first examples of a metal-to-insulator transition 
in which the insulating phase is generated by electron- 
electron correlations [SHE] . The Verwey transition also in- 
cludes sharp changes to the structure and magnetization 
[8]. Cooling through Ty, the unit cell structure changes 
from inverse-spinel to monoclinic and the magnetization 
exhibits a sharp decrease below Ty [8 . This coupling be- 
tween the electronic, magnetic, and structural properties 
makes magnetite very attractive for exploring tuning of 
magnetoelectric behavior in correlated materials which 
could lead to novel functionality. 

Previous demonstrations of electric field control of 
magnetic and metal-to- insulator transitions have resulted 
from various effects including current-induced breakdown 
of the insulating state [2j [3] [20] , field-induced changes of 
carrier concentration [4] [21] , field-induced strain gener- 
ated by growth on piezoelectric substrate (i.e. composite 
system) [5, 22-24 . The effect reported here is distinct 



from these previous categories. Unlike current-induced 
breakdown which is a highly non-equilibrium process, the 
present effect produces a true change in the equilibrium 
phase transition. We also find that this effect is not due 
to changes in carrier concentration, as shown by a sym- 
metric dependence of Ty on gate voltage. Finally, this 
effect does not rely on external strain provided by adja- 
cent layers. Thus, the electric field control of the Verwey 
transition represents a new type of electric field control 
in a highly correlated material. 

Fe304 films of 50 nm thickness are grown on double- 
side polished MgO(OOl) substrates using reactive molec- 
ular beam epitaxy (MBE) in ultrahigh vacuum (UHV) 
with a base pressure of 1 x 10 -10 torr. MgO substrates 
are first rinsed with de-ionized (DI) water. After load- 
ing into the MBE chamber, substrates are annealed at 
600° C for 45 minutes. A 10 nm MgO buffer layer is 
grown at 350° C via electron beam (e-beam) deposition 
from an MgO source [25]. Next, the Fes04 layer is grown 
at 200° C by depositing elemental Fe in a molecular oxy- 
gen partial pressure of 1.2 x 10 -7 torr. The Fe is evap- 
orated from a thermal effusion cell at a rate of ~0.13 
nm/min (for pure Fe). The single-crystal structure is 
verified through in situ reflection high energy electron 
diffraction (RHEED) and low energy electron diffraction 
(LEED), as shown in Figs. [T^l, [TJd, and^ inset. 0-20 
high resolution x-ray diffraction (HRXRD) scans exhibit 
a Fe 3 O 4 (004) peak near the MgO (002) substrate peak 
(Fig. [TJ^). Kiessig interference fringes indicate atomically 
smooth interfaces and verify the film thickness. 

Electrical properties of Fes04 films are characterized 
using standard dc four-point probe measurements (Fig- 
ure [l]i inset). Resistance values are obtained from 
current- voltage (/- V) curves, which exhibit linear depen- 
dence (Figure le) above 70 K. The temperature depen- 
dence of resistance (Figure [l]i, blue) exhibits a metal-to- 
insulator transition with a substantially higher resistance 
below 117 K, indicating the Verwey transition. Temper- 
ature dependence curves are measured as a function of 
increasing temperature, with the temperature stabilized 
for 10 min before a measurement is taken. For each tern- 
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FIG. 1. Characterization of Fes04 thin films, (a) and (b) 
are RHEED patterns for 50 nm Fe 3 4 on MgO(OOl) along the 
[100] and [110] in-plane directions, respectively, (c) HRXRD 
— 20 scans measured around the location of the MgO(002) 
peak with Kiessig fringes. Inset: LEED pattern with incident 
energy of 160 eV. (d) Temperature dependence of resistance 
measured by four-point probe (blue) and magnetization mea- 
sured by MOKE (red). The vertical dashed line indicates the 
Verwey transition. Inset: geometry for the resistance and 
magnetization measurements, (e) I-V curve for FesCU chan- 
nel at 115 K. (f) MOKE hysteresis loop for Fe 3 4 at 113 K. 
The right axis shows absolute magnetization based on SQUID 
measurements on corresponding samples. 

perature, measurements are repeated to ensure the tem- 
perature is stable. 

The magneto-optic Kerr effect (MOKE), with laser 
beam incident through the transparent MgO substrate, is 
used to characterize the magnetic properties of the Fe304 
films (812 nm wavelength, p-polarized, 45° angle of in- 
cidence). Figure If shows a typical longitudinal MOKE 
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FIG. 2. Electrical gating of Fes04 and manipulation of 
the Verwey transition, (a) A schematic of the sample de- 
vice structure. The dielectric layer consists of PMMA(900 
nm)/Al 2 03(50 nm)/MgO(10 nm). Positive electric field cor- 
responds to the application of positive voltage to the top gate 
electrode, (b) Temperature dependence of resistance for ap- 
plied electric fields of +1.8 MV/cm (blue), +0.9 MV/cm (or- 
ange), MV/cm (black), -0.9 MV/cm (green), -1.8 MV/cm 
(red). The arrows show TV for each electric field, which is 
summarized in the inset. 



hysteresis loop that exhibits large remanence and sharp 
magnetization reversal. The right hand axis of Fig. If 
displays the corresponding magnitude of the magnetiza- 
tion based on superconducting quantum interference de- 
vice (SQUID) magnetometry. The temperature depen- 
dence of the MOKE signal (Figure Id, red) exhibits a 
suppression of magnetization for temperatures below Ty . 
This behavior is characteristic of the Verwey transition. 

To apply electric fields to the Fes04 film, an insulating 
layer (PMMA/Al 2 3 /MgO) is deposited on top of the 
Fes04, followed by a metallic electrostatic gate (Pd/Ti), 
as shown schematically in Figure 2a. The devices are fab- 
ricated through several steps of evaporation using shadow 
masks. For the Fes04 layer, a narrow channel is produced 
with a width of 210 /im, creating a small active area to 
reduce the occurrence of pinholes and gate leakage. The 
Fes04 channel length is 4.2 mm and the gate length is 
3.3 mm. Alternate samples with Fe304 films covering a 
large area of the substrate produce similar results [26] . 
Pd(100 nm)/Ti(15 nm) contacts (for four-point probe) 
are e-beam evaporated through a shadow mask in a sep- 
arate system. Then a 10 nm MgO layer is grown on the 
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FIG. 3. Electrostatic gate dependence of resistance, (a) Gate 
dependent resistance for Fes04 at a temperature of 115 K. (b) 
and (c) show the symmetric and anti-symmetric components 
of the gate dependent resistance change, respectively. 



Fe 3 4 followed by a 50 nm A1 2 3 layer. PMMA is then 
spin coated onto the sample at 3000 rpm and cured un- 
der a vacuum environment at 170°C. The spin coating 
and baking sequence is repeated three times giving a fi- 
nal PMMA layer thickness of 900 nm. Finally, a shadow 
mask is used to grow the Pd(100 nm)/Ti(15 nm) top gate 
electrode. Typical gate leakage is 0.5 nA for electric fields 
of ±1.8 MV/cm. 

An electric field is produced by applying a voltage {Vq) 
between the gate electrode and the FesC^ film. Figure [2]o 
shows the temperature dependence of resistance for ap- 
plied electric fields of ±1.8 MV/cm (V G = ±60 V, blue), 
±0.9 MV/cm (V G = ±30 V, orange), MV/cm (V G = 
V black), -0.9 MV/cm (V G = -30 V, green), -1.8 MV/cm 
(Vg = -60 V, red) with corresponding colored arrows in- 
dicating Ty. The data clearly show that Ty varies as a 
function electric field, as summarized in the inset of Fig. 
[2J3. At zero electric field, Ty is 117 K. Strikingly, both 
positive and negative electric fields cause Ty to increase, 
indicating that the shift in Ty depends primarily on the 
magnitude of electric field as opposed to its sign. The 
maximum effect is observed for -1.8 MV/cm, where Ty 
increases to 119 K, giving ATy = ±2 K. At a similar 
electric field, the largest ATy we observe in our study is 
ATy = ±6 K for a large area sample [25] . The increase 
of Ty cannot be due to Joule heating because a heating 
artifact would appear as a reduction of Ty . We also rule 
out effects of irreversible sample change by measuring 
the zero electric field temperature dependence of resis- 
tance before and after taking the data in Fig. [2}d and 
no irreversible changes were observed. Finally, we ob- 
serve that temperature dependent resistance above Ty 
does not change with applied electric field, which shows 
that the metallic phase is insensitive to electric field. 

To gain further insight into the electric field effect, we 
perform a detailed study of the gate dependent resistance 
under isothermal conditions. Figure [3^i shows the resis- 
tance at 115 K as the electric field is swept between ±1.8 



MV/cm and -1.8 MV/cm. Consistent with the shift in 
Ty (inset Fig. J2Jd), the resistance increases for both pos- 
itive and negative electric fields and the effect is slightly 
larger for negative electric fields. To quantify the sym- 
metry of the electric field effect, we separate the change 
in resistance AR(E) = R(E) — R(0) into a symmetric 
part AR S (E) = [AR(E) ± AR{-E)]/2 (Fig. 3b) and 
anti-symmetric part AR A (E) = [AR(E) - AR(-E))/2 
(Fig. 3c). Comparing Fig. [3]o and [3]^, the symmetric 
part is up to 11 times larger than the anti-symmetric 
part. Because the change in carrier concentration is pro- 
portional to E (i.e. anti-symmetric), the small contribu- 
tion of ARa indicates that electric field control of the 
Verwey transition is not driven by a carrier concentra- 
tion effect. Instead, a symmetric effect can be driven 
by other interactions with the electric field. The pres- 
ence of an electric field will induce electric polarization 
given by P = XeE = — 1)^0^, where Xe is the electric 
susceptibility, eo is the permittivity of free space, and 
n is the relative dielectric constant of Fe304. The in- 
duced polarization will produce an energy contribution 
U = —\{PE) = —^(tv— l)eoE 2 that is symmetric in 
E. In addition, as Fe304 undergoes the Verwey transi- 
tion, the dielectric constant changes sharply with n being 
larger for the insulating state than for the metallic state 
(^ins > ^metai) El] • Thus, energy is lower for the insu- 
lating state than for the metallic state, which stabilizes 
the low temperature insulating state and causes Ty to 
increase. Therefore, this provides a macroscopic explana- 
tion for an electric field effect that is symmetric in E and 
produces an increase in Ty, consistent with experimental 
results. Further theoretical work is needed, including a 
microscopic model that can provide an explanation for 
the magnitude of the effect. In addition, a contribution 
from electric field induced strain could generate this sym- 
metry and should also be investigated [T8] . 

Since the Verwey transition in FeaC^ is a correlated 
phase transition that couples both the charge and mag- 
netic properties, it should be possible to tune magnetic 
properties with applied electric field. Figure [4^i shows 
MOKE hysteresis loops measured at 113 K with applied 
electric field of MV/cm (black) and -1.8 MV/cm (red). 
The absolute magnetization is determined by SQUID 
measurements (right axis of Figure |4^i) . An electric field 
of -1.8 MV/cm causes a decrease in the saturation mag- 
netization of 18%. Figure [4)3 displays the saturation 
magnetization as the electric field is swept between ±1.8 
MV/cm and -1.8 MV/cm. With the application of either 
positive or negative field, the magnetization decreases 
with a slightly stronger effect for negative fields. The 
magnetoelectric behavior is generated because the mag- 
netization has strong temperature dependence below Ty 
(Figure[l]i). When electric field is applied, the increase of 
Ty causes magnetization M to decrease because dM/dT 
is positive at T = 113 K; the intuitive picture is that the 
M vs T curve of Fig. [T]i shifts in temperature as Ty 
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FIG. 4. Electrostatic gate dependence of magnetization, (a) 
MOKE loops measured at 113 K with applied electric fields of 
MV/cm (black) and -1.8 MV/cm (red), showing a decrease 
in magnetization with the application of an electric field, (b) 
Magnetization as a function of electric field, demonstrating a 
magnetoelectric effect induced by electric field control of the 
Verwey transition. 



increases. Thus, the decrease of magnetization for both 
positive and negative fields (with slightly stronger effect 
for negative fields) is consistent with the electric field 
dependence of Ty (Fig. [2]o inset) and resistance (Fig. 
[3J. The change in magnetization as a function of elec- 
tric field is quantified by a magnetoelectric coefficient, 
&me = \ AM/AE\, where AM is the change in magneti- 
zation and AE is the change in electric field. Comparing 
the values at E = MV/cm and E = -1.8 MV/cm yields 
a value of olme — 585 ± 39 pT m/V. This is about an 
order of magnitude larger than in other single-phase mag- 
netoelectric materials, but not as large as some composite 
systems [26]. These results are compelling because they 
demonstrate a new method for generating magnetoelec- 
tric effects by controlling a correlated phase transition. 

In conclusion, we have demonstrated the electric field 
control of the Verwey transition in Fe3C>4 thin films. An 
electric field stabilizes the charge-ordered insulating state 
causing the Verwey transition temperature to increase. 
By manipulating a correlated phase transition that com- 
bines both charge and magnetic transitions, we realize a 
large and novel magnetoelectric effect. 
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